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Abstract 

The Millimeter Sky Transparency Imager (MiSTI) is a small millimeter- wave scanning telescope with 
a 25-cm diameter dish operating at 183 GHz. MiSTI is installed at Atacama, Chile, and it measures 
emission from atmospheric water vapor and its fluctuations to estimate atmospheric absorption in the 
millimeter to submillimeter. MiSTI observes the water vapor distribution at a spatial resolution of 0.5°, 
(— > , and it is sensitive enough to detect an excess path length of < 0.05 mm for an integration time of 1 s. 

By comparing the MiSTI measurements with those by a 220 GHz tipper, we validate that the 183 GHz 
measurements of MiSTI are correct, down to the level of any residual systematic errors in the 220 GHz 
measurements. Since 2008, MiSTI has provided real-time (every 1 hr) monitoring of the all-sky opacity 
distribution and atmospheric transmission curves in the (sub) millimeter through the internet, allowing to 
know the (sub)millimeter sky conditions at Atacama. 
' Key words: atmospheric effects — instrumentation: miscellaneous — site testing — submillimeter 
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o 

1. Introduction fore, necessary to make best use of a telescope under all 

conditions, and real-time monitoring of the atmospheric 
The major obstacle that leaves the millimeter (mm) condition is indispensable to execute it. Tippers operat- 
to sub-millimeter (submm) windows still uncxploitcd in ing at 220-225 GHz arc widely used at submm observato- 
ground-based astronomy is severe atmospheric absorption ries (e.g., CSO, ASTE, APEX, etc.) for sounding the sky 
. £h ' mainly caused by tropospheric water vapor. Moreover, opacity and assessing astronomical data quality. An all- 
small scale non-uniformity in its spatial distribution dis- sky mid-infrared (IR) imager to monitor clouds on the sky, 
$H ' torts wavefronts of astronomical signal, limiting the nat- similar to those commonly used at optical/near-IR obser- 
ural seeing at these wavelengths (e.g., Carilli & Holdaway vatories (Takato et al. 2003; Shamir et al. 2005; Suganuma 
1999). At the ALMA (Wootten & Thompson 2009) site, et al. 2007; Sebag et al. 2008; Miyata et al. 2008), is also 
precipitable water vapor (PWV) contents of the atmo- working at the ASTE site (4860 m in elevation; Pampa 
sphere have been the subject of several site testing cam- la Bola, Atacama, Chile; Ezawa et al. 2004; Ezawa et al. 
paigns using devices such as tipping radiometers, Fourier- 2008 for the ASTE facility). There has been, however, 
transform spectrometers (FTSs) and water vapor moni- no instrument to map the mm/submm emission of water 
tors (WVMs) operating at mm to submm wavelengths, vapor across the sky with high sensitivity. 
These experiments succeeded in modeling atmospheric It will be important to monitor not only the amount 
transparency at these wavelengths as a function of pre- of absorption by water vapor but also the amplitude of 
cipitable water vapor (Matsuo et al. 1998; Matsushita et small-scale (^10-100 m) fluctuations because they govern 
al. 1999; Paine 2000; Paine 2004; Pardo et al. 2001). the astronomical seeing at mm and submm wavelengths 
Another important fact from the viewpoint of telescope (see Altcnhoff et al. (1987) for single dish observations and 
operation at mm and submm wavelengths is that obscrva- Thompson, Moran, & Swcnson (2001) for intcrfcrometry) . 
tions above ~ 450 GHz usually require very good weather The turbulent layer of water vapor is known to be sal km 
conditions that are even at an excellent site like the above the ground level of the Atacama plateau (Delgado 
Chajnantor plain in Chile. Dynamic scheduling is, there- et al. 2001; Robson et al. 2002; Beaupuits et al. 2005). 
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The Millimeter Sky Transparency Imager (MiSTI) is a 
25-cm diameter telescope operating at 183 GHz installed 
at the ASTE site. It is designed to measure the emission 
from the atmospheric water vapor and its time variation in 
any arbitrary azimuth (Az) and elevation (El) direction. 
A picture and schematic drawing of MiSTI are shown in 
figure 1. and monitoring, and using it to verify the quality 
of ALMA data (i.e., data quality assurance). We have de- 
veloped MiSTI at the National Astronomical Observatory 
of Japan (NAOJ), and have installed and tested the in- 
strument at the high-altitude ASTE site. The 'first light' 
at the ASTE site was received in December 2007 (fig- 
ure 2), and stand-alone automatic/remote operations were 
started in March 2008. Since May 2008, a dedicated web 
site for MiSTI 1 has been opened to provide the real-time 
183 GHz sky images to radio astronomers who already 
have (and will have) telescopes in Atacama. 

In this paper, we present the instrumental overview of 
MiSTI and its initial results; more detailed analyses of the 
sky measurements will be presented in a subsequent paper. 
Requirements on the instrument and methods for measur- 
ing water vapor are introduced in Section 2. The system 
overview of MiSTI is described in Section 3. Application 
of the MiSTI as an all-sky opacity monitor is presented in 
Section 4. Sections 5 and 6 devote discussions and sum- 
mary, respectively. 

2. Requirements and Methodology 

2.1. Science Goals and Requirements 

One of the scientific goals of MiSTI is to map and moni- 
tor the mm emission of water vapor across the sky to use it 
for dynamic scheduling of astronomical observations and 
data quality assurance. This allows an observer at submm 
observatories to choose a transparent region of the sky to 
be observed and to make the observing schedule more ef- 
ficient. This requires high sensitivity to water vapor and 
the ability to rapidly point to arbitrary positions of the 
sky, since all-sky mapping with a single beam is generally 
time-consuming. 

The other goal is to understand small-scale (~10— 
100 m) structure of water vapor distribution ~sl km above 
the ground level. This can only be done with a radiometer 
by quickly measuring the amount of water vapor at ~«2- 
9 points within a small patch of the sky, a l°-5° square 
region for instance. The measurement is on an experi- 
mental basis and will be very challenging. But it is im- 
portant because the small-scale fluctuation of water vapor 
is believed to affect natural seeing of submm single-dish 
telescopes and interferometers. This requires a spatial res- 
olution of ~ 0.5° and a quick scanning or wobbling ability. 
Furthermore, sensitivity to water vapor and a high stabil- 
ity receiver are necessary to measure the tiny fluctuations 
of water vapor even for short integration times of ~1 s. 
Portability of the instrument is also required since exper- 
iments may be performed at different observatories and 
sites. 

1 http : //aste-www .mtk .nao .ac.jp/~misti/ opacity .html 



The required sensitivity to water vapor and diffraction- 
limited spatial resolution depend on observing frequency. 
In the next section we examine the optimal observing fre- 
quency. 

2.2. Selection of Observing Frequency 

One of the most effective ways to detect small vari- 
ations of water vapor is to observe a water line. In 
the wavebands which have a significant water line, wa- 
ter vapor emission increases the closer to the line centers. 
Furthermore, the amount of variation in brightness tem- 
perature for a unit variation in the amount of water va- 
por also increases, which improves the sensitivity to tiny 
changes in PWV. Water vapor emits several lines at cen- 
timeter to submm wavelengths (e.g., 22, 183, 325 GHz). 
We choose the 183 GHz line for the following reasons. 
The 22 GHz line is always optically thin in the dry con- 
ditions at Atacama, and hence the sensitivity to varia- 
tions of PWV is lower than at 183 GHz, as suggested by 
Wiedner et al. (2001). Moreover, receivers at 325 GHz 
generally have higher noise temperatures than 183 GHz. 

In order to optimize the radio frequency that is effi- 
cient to detect small variations of water vapor, we use an 
atmospheric transmission model ATM (Pardo et al. 2001) 
to characterize the sky brightness temperature and excess 
path length in the 183 GHz band. We assume an alti- 
tude and vertical pressure distribution similar to those in 
Atacama. PWV is set to 0.8, 1.4, and 2.8 mm, which 
are the first, second and third quartiles (i.e., 25%, 50% 
and 75%) of PWV in Atacama, respectively (Radford & 
Holdaway 1998). 

In the top panel of figure 3, we show the variation in 
brightness temperature per unit change in excess path 
length, ST/SI, as a function of frequency. The middle 
panel shows the uncertainties in sky temperature mea- 
surement, AT, as a function of frequency. The uncertainty 
is expressed as AT = (T RX + T° ky )/ ^/Bt intcg , where T* ky 
is the sky brightness temperature, and we assume a re- 
ceiver noise temperature Trx = 1000 K, the band width 
B = 2 GHz (DSB) and the integration time ii n t og = 1 s. If 
we consider the former and latter as the signal and noise, 
respectively, then effective sensitivity to water vapor emis- 
sion is proportional to a ratio between ST j SI and AT. The 
bottom panel of figure 3 shows the minimum detection 
limit to excess path length (EP). The local minima of the 
detection limit appear at ^ op t ~ 183.3±(3-4) GHz, and 
the position of the local maxima slightly moves outward 
from the water line center with increasing PWV. This is 
because water vapor emission at frequencies close to the 
line center is rapidly saturated with increasing PWV and 
therefore it becomes relatively insensitive to variation of 
PWV. We expect that the water line is almost saturated 
for PWV > 2.8 mm at 180-186 GHz. This situation is 
more likely to occur when the telescope points at low el- 
evation angles. We prefer the frequency band in which 
ST j SI is sufficiently large and the water vapor emission 
is not saturated over a broad range of PWV. We there- 
fore adopt 183. 3 ±4.5 GHz as the observing frequency of 
MiSTI. 
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We find that a 183 GHz receiver operating at room tem- 
peratures with Trx ~ 1000 K and B = 2 GHz is sensitive 
enough to detect an EP of < 0.05 mm for 1 s integration. 
Note that if a difference in EP at two different points of 
the sky separated by a distance of 10 m is 0.05 mm, a 
wavefront falling on a 10 m telescope aperture can tilt at 
~ 1" (see the right axis of the bottom panel in figure 3). 

3. The System 

3.1. Telescope 

MiSTI is a 1-m tall, portable mm-wave telescope with a 
25 cm primary mirror. The specification and system block 
diagram are shown in table 1 and figure 4, respectively. 
The optics and receiver of MiSTI is on an alt-azimuth 
mount. The telescope boresight can move in the ranges 
of Az = -175° to +175° and El = 0° to 95°. For the 
measurements of small-scale structure of water vapor, the 
elevation drive unit equipped with a main reflector, sec- 
ondary and tertiary mirrors can be wobbled at an angle of 
30' at 1 Hz. This allows to measure PWV at two points 
on the sky, «9 m apart from one another if the water va- 
por layer is 1 km above the ground. The elevation drive 
unit is outside the receiver cabin and it is covered with a 
0.5-mm-thick GORE-TEX membrane to avoid desert dust 
and snow (see figure 1). The instruments for measurement 
and control, except for a control computer, are inside the 
lower mount. 

MiSTI is installed on a frame firmly welded on the top 
of a container at the ASTE site. The latitude, longitude, 
and altitude are N22°58'18", W67°42'10", and 4860 m, 
respectively. The infrastructure (e.g., power supply and 
internet access) is provided by the ASTE facility. 

3.2. Optics 

The optics of MiSTI is a dual-offset modified-Nasmyth 
design, and consists of a main paraboloidal reflector (MR), 
a secondary hyperboloidal reflector (SR), two ellipsoidal 
mirrors (M3 and M4) and a feed horn. The main reflec- 
tor is an 250-mm offset paraboloidal mirror, and it has 
a focal length of / = 245.0 mm. The main reflector is 
under- illuminated: the edge taper level is T c = — 20 dB, 
which provides an angular resolution of 30 arcmin at 
183.3 GHz and low (^1%) spillover loss. The beam shape 
was checked by raster-scanning the Sun. The sky emis- 
sion collected by the primary, secondary and tertiary mir- 
rors is guided into the receiver cabin along the elevation 
axis, and it is fed into the mixer through the corrugated 
feed horn after being focused by the forth mirror. For 
antenna temperature calibration, MiSTI has two offset- 
par aboloidal mirrors, each of which moves and chops the 
ray at a beam waist between M3 and M4. Then, each 
chopping mirror bends the beam from the feed to couple 
it with radio absorbers at different temperatures. All of 
the optical elements are optimized to 183.3 GHz on the 
basis of Gaussian optics. 



3.3. Front-end and back-end system 

The front-end system consists of a diode mixer, a low- 
noise amplifier (LNA), a bandpass filter (BPF), and a sec- 
ond amplifier, all of which are operated at room tempera- 
ture (i.e., not cooled). We use a subharmonic mixer, WR- 
4.3SHM (Virginia Diode, Inc.), which has a DSB noise 
temperature and a conversion loss of 600 K and 5.5 dB, 
respectively, over the operating frequencies 175-220 GHz. 
The local oscillator (LO) signal is supplied from a Gunn 
oscillator (Wise Wave Technologies, Inc.). The Gunn oscil- 
lator has no phase-lock loop, since we found that the LO 
frequency always falls within a narrow frequency range 
(<100 MHz) centered at 183.3 GHz. Furthermore, we 
made sure that the uncertainty in opacity estimates is very 
small (< 1%) for any conditions at PWV = 0.2-2.0 mm 
even if the LO frequency shifts by 100 MHz from the wa- 
ter line center. This is because a shift of LO frequency can 
increase power received by one of the sidebands and de- 
crease power received by the other sideband, which almost 
compensates the total detected power. An attenuator in- 
serted between the oscillator and mixer is used to optimize 
the power of the LO signal. The intermediate frequency 
(IF) is centered at 4.5 GHz with a bandwidth of 1.0 GHz. 
The receiver noise temperature is Trx ~ 900 K. The biases 
of the amplifiers and oscillator are fed by a direct-current 
(DC) power supply (E3631/Agilent Technologies, Inc). 

We employ a power sensor and meter (E9300A and 
E4418B /Agilent Technologies, Inc.) as an IF total power 
detector and its readout instrument. The readout sig- 
nal detected with a power sensor is modulated (i.e., AC- 
coupled) to separate the signal from the low-frequency 
noise. This significantly improves suppression of l//-type 
noise. 

Figure 5 shows the Allan variance of the output of the 
power meter. We measured the output at 196 GHz for 
3600 s only turning on the devices used for receiver oper- 
ation and data acquisition. The feed was coupled with a 
radio absorber at a constant temperature of 297 K. The 
power meter data were normalized by the power averaged 
over the measurement. The Allan variance shows a ther- 
mal behavior down to ~100 s where it starts to flatten at 
a level of o^ llan ~2x 10 -9 . This indicates a gain stability 
of AG/G w V2 x 10" 9 - 5 x 10~ 5 for a time duration of 
< 100 s. When the devices for full operation (including 
motor drivers, hot load system, and fan) are powered on, 
the gain stability slightly degrades probably due to elec- 
tromagnetic interference from those devices. We, however, 
find (Taiian ~ 3 x 10 -9 for 50 s, suggesting that the gain is 
still quite stable (AG/G < 10~ 4 ) for 50 s. 

3.4. Intensity Calibration 

For measurements of antenna temperature of water va- 
por emission, the dual-load calibration method is ap- 
propriate for intensity calibration since we do not need 
prior information about atmospheric temperature, which 
is often difficult to know. MiSTI is equipped with dual- 
temperature loads in its receiver cabin. One is op- 
erating at an ambient temperature, and the other is 
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heated. Temperatures of both of the loads are moni- 
tored. The beam radiated from the feed is folded by 
rotating paraboloidal mirrors, which focus the beam on 
temperature standards. By making the radiated area on 
the loads as small as possible (beam radius on the loads, 
w = 3.6 mm), the load temperature measurements are eas- 
ier and more accurate. 

We use Radio Absorbing Material (RAM; Thomas 
Keating Inc.) for temperature standards. We insert 
a high-precision temperature sensor module (LM35CA; 
National Semiconductor) in each of the RAMs to monitor 
their temperatures. The absolute accuracy is 0.4°C over 
an operating temperature range from —40 to +110°C. 

We have paid special attention to keeping the temper- 
ature of the hot load as stable as possible. The hot load 
system consists of a RAM, a flat heater (SAMICON 230; 
Sakaguchi E.H Voc Co.), a temperature control module 
(TTM-000; Toho Electronics Inc.) and a small circuit for 
safety and temperature control. The radio absorbers are 
settled on an aluminum holder, which has high thermal 
conductivity and is used for heating the RAM from its 
back and sides to make the surface temperature of the ab- 
sorber uniform. We filled the gap between the absorber 
and the holder with silicone grease to improve the thermal 
conductivity. The window to the RAM is covered with a 
foam polystyrene sheet (1 mm in depth) to keep the sur- 
face temperature of the RAM high and constant. The 
temperature controller monitors the RAM temperature 
with a thermocouple, and it keeps the temperature con- 
stant with a proportional-intcgral-derivative (PID) con- 
trol algorithm. We tested the hot load system and found 
that the system can stably and safely heat up the absorber 
up to «110°C. 

An optional cold temperature standard can be available. 
The receiver cabin has a space in which a shallow dewar 
flask with RAM at the bottom is set. The flask can be 
filled with 1 litter of liquid nitrogen to cool the RAM. 
The chopper mirror can bend the beam and illuminate 
the cooled RAM. The boiling point of liquid nitrogen at 
4860 m altitude (air pressure of 570-575 hPa) is ss 73 K. 

For measurements of atmospheric transparency by sam- 
pling powers from the sky at several or more elevation 
angles, commonly referred to skydip measurements, the 
one-load calibration method is still useful. In the one-load 
method, an atmospheric temperature is assumed to be 
equal to the ambient load temperature. The formulation is 
given in Section 4.1. Note that if an atmospheric tempera- 
ture is lower than an ambient temperature by 10 K, MiSTI 
may underestimate the opacity by ~10% for a opacity of 
~0.4, which is approximately a median optical depth of 
the atmosphere in Atacama at 183±4.5 GHz. 

3.5. Control System 

The hardware for control and measurement (e.g., motor 
drivers, power meter, digital multimeter) was put inside 
the mount box of MiSTI as mentioned above, whereas a 
dedicated personal computer (PC) is placed in the con- 
tainer right below the telescope. The control software is 
installed on the PC, and it controls the azimuth and el- 



evation motors, the chopping mirrors, the power meter 
for IF power measurements, and a digital multimeter for 
temperature measurements via a local area network and a 
GPIB interface (GPIB-ENET/100, National Instruments 
Inc.). See also figure 4. The control software runs on the 
Linux operating system, and is written in C and Perl lan- 
guages. Observations are performed by running a script 
that has azimuth and elevation coordinates for each inte- 
gration and is prepared in advance. The software reads 
the script and sends to the motor drivers the number 
of pulses that corresponds to the pointing coordinates. 
When the telescope points to the intended direction of 
the sky, the software starts to acquire and save the out- 
put from the power meter. Simultaneously the software 
measures and records the load temperatures. Currently, 
observing scripts for standard scans (e.g., skydips, raster- 
scans and all-sky maps) and the associated reduction soft- 
ware are available. 

4. The Sky Opacity Monitor 

MiSTI is operating as the ASTE Sky Monitor, which 
is open to the public through a web site 2 . The monitor 
provides useful information for constructing (sub)mm ob- 
serving strategies at the Pampa la Bola and Chajnantor 
sites. The ASTE Sky Monitor includes two monitors: (1) 
the all-sky opacity monitor, which shows the 183-GHz 
opacity distribution over a whole sky and (2) the atmo- 
spheric window monitor, which offers transmission curves 
of the atmospheric windows. Here we provide description 
of those monitors. 

4-1- All-Sky Opacity Mapping 

The radiometer measures the power from the sky at 
each azimuth and elevation, P s \ l y(Az, El). The integra- 
tion time at each Az and El is 0.4 sec, which provides a 
receiver noise level of ~0.05 K (rms). For intensity cali- 
bration, we employ the one-load chopper wheel method at 
the moment and will update it to the dual-load method. 
The radiometer measures the power from a load, P ca \(t), 
once every 50 integrations corresponding to a time dura- 
tion of ~30 sec. The load is set inside the receiver cabin 
and its temperature is not actively controlled but moni- 
tored. The all-sky mapping typically takes ~40 min, and 
the monitor is updated hourly. 

When the absorbing layer of the sky is approximately 
expressed as a plane parallel slab with a finite thickness 
in height and the atmosphere and load temperatures can 
be regarded to be the same, the relationship among P s ky, 
P ca i, Az and El can be described as 



In 



1 



P sky (Az,El;t) 



Pcal(i) 



-tq sec Z + b, 



(1) 



where Z is the zenith angle defined as (90 — El) deg, and 
To and b are positive constants. The constant To is called 
zenith opacity, which is the optical depth of the sky ab- 
sorbing layer at the zenith (El = 90°). We derive the 
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constants tq and b by fitting the data to a linear function 
of secZ. The opacity in an arbitrary direction (Az,El) is 
then described as 



r(Az,El) = b-ln 



1 



P s]iy (Az,El;t) 



(2) 



actually measured with the tipper. Again, the correla- 
tion is very good (the correlation coefficient is 0.958). 
When the opacities are relatively low (tq 20 < 0.08, or 
70% of the data points in our comparison), the transfcrcd 
220 GHz opacities are consistent with the tipper's opac- 



ities (Tq 



220 



r 220 MiSTI 
r 



). This suggests that the 183 GHz 



(3) 



We define the opacity fluctuation 5t(Az,EI) as 

dr(Az,El) = T{MEl) -l 0SecZ . 

To sec Z 

An example of the all-sky measurements and results is 
presented in figure 6. The scanning pathway is shown in 
figure 7. 

4-2. Atmospheric Modeling 

The 183 GHz opacity (t 183 ) is a good measure of PWV 
when the measured sky opacity < 1 as discussed in § 2.2. 
MiSTI exploits the advantage of the 183 GHz spectrum, 
and provides not only PWV but also the transmission 
curves of the atmospheric windows by introducing an at- 
mospheric transmission model, AM (Paine 2004). AM is a 
program for producing model atmospheric emission, ab- 
sorption, and phase delay at frequencies between a few 
GHz and a few THz. The model that we use here is 
optimized for the Chajnantor site (am configuration file, 
Chajnantor_lowres . amc). The model does not consider 
liquid phase of water. The dominant parameters that con- 
strain the atmospheric spectra (or opacity at an arbitrary 
frequency) are PWV and ground temperature, T gr , which 
allows us to estimate PWV from measurements of t 183 effect from water droplets contained in visible clouds, 
and T gr . We use the atmospheric model to calculate Tq 83 
for each PWV and T gr , and then we search for the PWV 
that accounts for the measured To and T gr using the least- 
X 2 method. Once PWV is determined, we can then es- 
timate the opacity spectrum over the mm/submm bands 
using AM. 



measurements of MiSTI are correct, down to the level of 
any residual systematic errors in the measurements of the 
220 GHz tipper. 

In contrast, when the opacities arc relatively high 
(tq 20 > 0.1), the actual 220 GHz measurements clearly 
exceed the transfered opacities. The best-fit line to all of 
the data is expressed as r 220 = 1.1432 r 220MiSTI + 0.002, 
suggesting a ^14% excess of Tq 20 to the MiSTI measure- 
ments on average. This might be due to a calibration error 
in load-sky measurements at 183 GHz, which decreases 
the 183 GHz optical depths, as suggested in Section 3.4. 
However, the calibration scheme, the load material and 
temperature used in MiSTI are the same as those used in 
the tipper. It is thus unlikely that the calibration method 
itself largely affects the discrepancy. 

Another more likely origin of the excess can be the effect 
of water droplets on the optical depth measured with the 
220 GHz tipper. It is known that not only water vapor 
but also liquid water contribute to the optical depth at 
220 GHz whereas the 183 GHz opacity depends almost 
only on pure water vapor (e.g., Matsushita et al. 2000). 
It is thus likely that the liquid phase of water more affects 
on t, 220 than Tn 83 . In the next section, we discuss the 



5. Discussions 

The sky data obtained with the IR cloud monitor at 
the ASTE site can be used to examine the effect of the 
liquid component of water on the correlation between 



4- 3. Comparing Measurements at 220 GHz 

To assess how correctly the opacity measurement and 
modeling (i.e., frequency transfer of atmospheric opac- 
ity estimates) work, wc compared the zenith opacities 
measured at 183 GHz with those at 220 GHz. The 
220 GHz data were retrieved from the data archive of 
the tipper at the ASTE site (Kohno et al. 1995), which 
have been cross-checked by side-by-side experiments in- 
volving the 220 GHz tipper and other instruments such 
as the NRAO 225 GHz tipper (Radford & Holdaway 
1998) and an FTS (Matsuo et al. 1998; Matsushita et al. 
1999; Matsushita et al. 2000). We estimated the 220 GHz 
opacities from our 183 GHz measurements by modeling 
the atmospheric transmission, and compared them with 
the actual 220 GHz measurements. 

Figure 8a shows the opacity correlation between 183 
and 220 GHz. The data were collected during March to 
August and October to December in 2008. The correla- 
tion is excellent and the correlation coefficient is 0.957. 
The best-fit linear function is r 220 = 0.186t 183 + 0.011. In 
figure 8b, we compare the 220 GHz opacities which are 
estimated from the 183 GHz measurements, with those 



-183 



and 



r 220 



The cloud monitor consists of a gold- 



plated Cassegrain-type mirror, a commercial mid-IR cam- 
era which has 320 x 240 image pixels, and a Linux based 
computer for control and data gathering. The cloud mon- 
itor can automatically take and store a fish-eye image of 
the sky every 5 minutes (figure 9). The camera is sensitive 
to the wavelength range of 8-12 /im. Before installing the 
camera to the ASTE site, we measured its characteristics 
in the laboratory by taking images of an objects with vari- 
ous temperature (T bj), and we confirmed that the output 
of each pixel is quite proportional to 

/>12/^m 

/ B x {T ohi )d\ (4) 

in the range of T Q bj = 255-330 K, where B\(T) is the 
Planck function. In each frame of the sky image, the IR 
emission from a black-painted aluminum block placed at 
the edge of the main mirror is also taken simultaneously 
(see figure 9). The temperature of the block is always 
monitored, and its emission is used as the brightness stan- 
dard when the emissivity distribution of the sky is calcu- 
lated. The obtained sky emissivity is always above 0; Even 
when there are no clouds, the emissivity is uniform and its 



G 



Tamura et al. 



typical value is 0.23-0.33, or 0.25-0.4 in optical depth at 
these IR wavelengths (tir). This is mainly because there 
is an emission band of atmospheric ozone in this wave- 
length range. The measured accuracy is also limited by 
the accumulation of dirt on the mirror surface. More de- 
tailed description about the system can be found in the 
papers reported by Takato et al. (2003), Suganuma et al. 
(2007), and Miyata et al. (2008). 

Figures 10a and 10b show the comparison of optical 
depths at mm wavelengths and that at mid-IR during the 
period from April to August in 2008. These figures contain 
data that were taken only when all three environmental 
monitors at the ASTE site (i.e., MiSTI, tipper, and IR 
cloud monitor) were working. The time resolution of the 
data in figure 10 is smoothed to be 1 hour, which corre- 
sponds to 1.5 times the time resolution of MiSTI. When 
t 220 < 0.05 or t 183 < 0.4, r IR = 0.25-0.4, indicating that the 
sky was clear, while when the optical depths at mm wave- 
lengths are large (t, 220 > 0.05 or Tq 83 > 0.4), tir is large 
but also shows large scatter (figures 10a and 10b). Tq 83 
changes by an order of magnitude even under clear sky 
conditions, suggesting that the observed frequency used 
by MiSTI can give us critical information to select the best 
conditions for submm observations. Figure 10c presents 
comparisons between Tq 20 and Tg 83 under clear sky (i.e. 
tir < 0.4, dots) and cloudy conditions (tir > 0.4, open 
circles). These comparisons reveal the followings: (i) Tq 20 
is always greater than 0.05 when the sky is cloudy, (ii) 
the correlation is much tighter under clear sky conditions, 
and (iii) Tq 20 under cloudy sky shows large scatter and 
is even larger than the value expected from the correla- 
tion under clear sky. These results imply that the liquid 
component of water affects Tq 20 more than Tq 83 , as dis- 
cussed by Matsushita et al. (2000). The correlation under 
clear sky shown in figure 10c shows a bend at (tq 83 ,Tq 20 ) 
= (0.4, 0.06) and the slope gets steeper in the regions of 
larger Tq 83 , which is also seen in figure 8a. This may also 
be explained by more significant contribution of liquid wa- 
ter to t 220 . 

6. Summary 

We have developed a small mm- wave telescope equipped 
with a 25-cm diameter dish, named the Millimeter Sky 
Transparency Imager (MiSTI) . It measures emission from 
atmospheric water vapor at frequencies of 183.3 ±4.5 GHz, 
which maximize the sensitivity to tiny changes in the 
PWV content of the atmosphere. MiSTI observes the wa- 
ter vapor distribution at a spatial resolution of 0.5°, and 
it is sensitive enough to detect an excess path length of 
< 0.05 mm for an integration time of 1 s. By compar- 
ing the measurements of MiSTI at 183 GHz with those 
from the 220 GHz tipper, we validate that the 183 GHz 
measurements of MiSTI are correct, down to the level 
of any residual systematic errors in the measurements at 
220 GHz. When the mm opacities are large (tq 83 > 0.4), 
optical depths measured at 220 GHz start to exceed those 
at 183 GHz. Given the comparison between the mm 
and IR opacities, this excess is likely attributed to wa- 



ter droplets in the sky, which contribute more largely to 
220 GHz opacity than 183 GHz. Currently, MiSTI pro- 
vides real-time (every 1 hr) monitoring of the all-sky opac- 
ity distribution and atmospheric transmission curves in 
the (sub)mm through the internet, allowing us to know 
the (sub)mm sky conditions at Atacama. It is important 
to investigate ^10-100 m scale structure of water vapor 
distribution because the small-scale fluctuations of the wa- 
ter vapor affect the natural seeing of submm single-dish 
telescopes and interferometers. We plan to update the 
functionality of MiSTI to further increase its sensitivity 
and measurement efficiency of such fluctuations. 
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Table 1. The specifications of MiSTI. 



Parameter Value 



Main Reflector 



Optics 




offset paraboloid 


Diameter 


D 


250 mm 


Confocal distance 




-30 m 


Beam size (HPBW) 


A6 


30 arcmin 


Surface accuracy 




< 20 fj,m 


Edge taper 


T e 


-20 dB 


Predicted spillover efficiency 




0.99 


Receiver 


Receiver system 




heterodyne radiometer 


Observing frequency 


Vobs 


183.3 ±4.5 GHz 


Intermediate frequency 


I/IF 


4.0-5.0 GHz 


DSB band width 


B 


2.0 GHz 


Receiver noise temperature 


Trx(DSB) 


w900 K 


Gain stability 


AG/G 


<io~ 4 


Misc. 


Calibration method 




one/dual- load chopper wheel method 


Mount 




alt-azimuth mount 


Rotation speed 




10 deg s" 1 


Beam switch 




wobbling main reflector (1 Hz, 30') 





Fig. 1. {left) A picture of the Millimeter Sky Transparent Imager (MiSTI) installed at the ASTE site. The elevation drive unit 
equipped with the primary, secondary, and tertiary mirrors is covered with a thin GORE-TEX membrane to avoid desert dust and 
snow, (right) Schematic drawing of MiSTI. The membrane and its support structures are not shown. The components in cyan are 
optical mirrors, and the line in red shows the optical axis. The positions of the azimuth and elevation axes are noted by "AZ" and 
"EL", respectively. 
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l^ 1 









Fig. 2. {left) The first-light image of MiSTI. The target object is the ASTE 10-m dish. The receiver of MiSTI was tuned at a 
local-oscillator frequency of 196 GHz, where the atmospheric transmission is higher than at 183 GHz. The brighter colors represent 
higher intensities at 196 GHz. We also show the optical image of ASTE {right). 





180 



200 



Frequency (GHz) 



Fig. 3. {top) The change in brightness temperature of the sky per unit change in excess path (EP), ST/81. The dashed, solid, 
and dotted curves show ST/81 for the first, second and third quartiles (i.e., 25%, 50% (median), and 75%) of PWV at Chajnantor, 
respectively. The hatched regions indicate the frequencies at which MiSTI observes the sky. {middle) The uncertainty of sky 
brightness measurement as a function of frequency when assuming the integration time ti ntcg = 1.0 sec, the receiver noise temperature 
Trx = 1000 K and the double-sideband bandwidth B = 2 GHz. {bottom) The effective sensitivity (5<r) to EP length I as a function 
of frequency. The scale of the corresponding detection limit of wavefront tilt (pointing jitter for a 10-m telescope) is indicated in the 
right axis. See Section 2.2 for details. 
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Fig. 4. (left) The system block diagram of MiSTI. The boxes outlined by dashed lines indicate the mount and instrument room 
(grey box) , the azimuth unit (pale- green) , the elevation unit (yellow) and the chopping mirrors (orange) , and they are independently 
driven by servo motors and their drivers (MDs). Emission from the sky is collected by main reflector (MR), sub-reflector (SR) and 
the third and forth mirrors (M3, M4). Then, the emission is coupled by a feed horn, and detected by the receiver and the power 
sensor. The reference signal at 183.3 GHz is generated from a Gunn oscillator, and coupled at the mixer. For flux calibration, the 
chopping mirrors (CHP) bend the beam fed by the receiver, and couples the beam with loads (TK-RAM), whose temperatures are 
controlled and/or monitored by the temperature controller and digital multimeter (DMM). The power meter, the digital multimeter, 
and the motor drive units are controlled from a Linux machine via a local area network at the ASTE site, (right) A picture of 
the inside of MiSTI. The dashed lines in grey, pale-green, and yellow indicate the mount, the azimuth unit, and the elevation unit, 
respectively. The primary mirror in the picture is temporarily covered with a blue protective film. 




1 10 100 1000 

Time (s) 

Fig. 5. The Allan variance of receiver output at 196 GHz. The feed was coupled with a radio absorber at a temperature of 297 K. 
We measured the data with turning on only devices for receiver operation and data acquisition (i.e., power meter, DC supply, and 
network devices; black curve) and all devices for full operation (i.e., power meter, DC supply, network devices, motor drivers, hot 
load, and fan; grey curve). 
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Fig. 6. An example of the all-sky opacity image taken with MiSTI on 2010 July 1. (a) The total power which the receiver measured 
on the sky (black dots) and on the room load (red dots), as a function of time in UTC. (b) The fit. The best-fit linear function is 
shown in the red line. The zenith opacity at 183 GHz is Tq 83 = 0.066, which is almost the best conditions at Pampa la Bola. (c) 
Calibrated time-stream, (d) All-sky opacity map at 183 GHz in the fish-eye lens view defined as Equation 2. The sky hemisphere is 
orthographically projected to the plane of the paper. The north is top and the east is left. The solid and dotted circles indicate the 
horizon (El = 0°), El = 20,40,60, and 80°. (e) The fluctuations of the all-sky opacity map, 8t(Az,EI), which is define as Equation 3. 
In this measurement, the opacity typically fluctuates by 8% (rms) at El > 60°. See the text for details. 



z (zenith) 




Fig. 7. A schematic picture of the scanning pathway of the MiSTI's pointing used in the all-sky opacity measurement. The sphere 
represents the sky in the horizontal coordinate system, and the white disk perpendicular to the z axis represents the horizon. The 
open circles on the sphere indicate the direction which MiSTI points to. The positions on the sky hemisphere are selected, such that 
the corresponding positions orthographically projected to the horizon plane (x,y) are on regular grid points (xi,yj), where i and j 
arc integer. The azimuthal and zenith angles are then given by Az = tan -1 (yj/xi) and Z = sin -1 yjx^ + y?. The thick and thin 
arrows connecting the open circles show the scanning path on the sky and on the horizontal plane (or an image shown in figure 6d). 
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Fig. 8. (a) The comparison of zenith opacities measured with the 220 GHz tipper (tq 20 ) and MiSTI (tq 83 ). The error bars on 
each plot are calculated from rms of tipper outputs during each of the MiSTI measurements («41 min). The solid line shows 
the best-fit linear function, t 220 = 0.186r 183 +0.011. (b) The comparison of zenith opacities measured with the 220 GHz tipper 
(tq 20 ) and those estimated from the MiSTI 183 GHz measurements ( r 220MiSTl^ The so i^ ) me shows the best-fit linear function, 
r 22 ° = 1. 132 T 220 MiSTI + 0.002. The dashed line indicates t 220 = r 220 MiSTI Sec § 4 3 for morc dctails 



Fig. 9. Examples of emissivity distribution derived from the images taken with the IR cloud monitor at the ASTE site, (a) under 
clear sky condition (ttr ft: 0.25), and (b) under cloudy weather (tir si 0.5). Two aluminum blocks "Rcfl" and "Ref2" are placed at 
the edge of the main mirror to be used as brightness standards, but only "Refl" was used in the analyses shown in this paper. Note 
that the shadows by these blocks, sub-reflector and its pillars are masked in these images. 
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Fig. 10. Comparisons among the data taken with the tipper working at 220 GHz, MiSTI and the IR cloud monitor during the 
period from April to August in 2008; only data that were taken simultaneously with these three monitors are included. The time 
resolution of the data is smoothed to be 1 hour, (a) Comparisons between the 220 GHz optical depth toward zenith (tq 20 ) with 
the tipper and the IR optical depth averaged over the regions where the elevation is greater than 60° (tir) derived from the images 
with the IR cloud monitor. There are 2279 data points in total. The dashed line indicates the border between clear and cloudy 
conditions (tir =0.4). (b) Comparisons between the 183 GHz optical depth toward zenith (tq 83 ) with MiSTI and the tir. The 
number of the data points is the same as (a), (c) Comparisons between Tq 83 and tS 20 when the sky is clear (tjj^ < 0.4, black dots) 
and cloudy (-TrR > 0.4, red open circles). There are 2108 and 171 data points for the tir < 0.4 and tir > 0.4 conditions, respectively. 



